Europhysics Letters 



PREPRINT 



> 



Event-based simulation of single-photon beam splitters 
and Mach-Zehnder interferometers(*) 



H. De Raedt^, K. De Raedt^ and K. Michielsen 



3 



^ Department of Applied Physics, Materials Science Centre, University of Groningen, 
- - - Nijenborgh 4, NL-9747 AG Groningen, The Netherlands 

'iO I ^ Department of Computer Science, University of Groningen, Blauwborgje 3, NL-9747 

■ AC Groningen, The Netherlands 

(N 
Ch . 

^ PACS. 02.70.-c - Computational Techniques. 

PACS. 03.65.-W - Quantum Mechanics. 

(N 



Abstract. - We demonstrate that networks of locally connected processing units with a prim- 
itive learning capability exhibit behavior that is usually only attributed to quantum systems. 
, We describe networks that simulate single-photon beam-splitter and Mach-Zehnder interferom- 

eter experiments on a causal, event-by-event basis and demonstrate that the simulation results 
are in excellent agreement with quantum theory. 
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^ il Introduction. - The objective of the research reported in this paper is to demonstrate 

that locally-connected networks of processing units with a primitive learning capability are 
^ sufficient to simulate, event-by-event, the single-photon beam splitter and Mach-Zehnder inter- 
j5 ' ferometer experiments of Grangier et al. [1] . This is one of the basic experiments in quantum 
physics [1,2] that has not been simulated in the event-by-event manner in which the experi- 
• • . mental observations are actually recorded [3] . Although quantum theory gives us a recipe to 
. 5^ ' compute the frequency for observing different types of events, it does not describe individual 
. events [2,4,5]. Reconciling the mathematical formalism (that does not describe single events) 
?H ' with the experimental fact that each observation yields a definite outcome is often referred 
. 5^ 1 to as the quantum measurement paradox. This is a central, fundamental problem in the 
foundation of quantum theory [5-7] . 

From a computational viewpoint, quantum theory provides us with a set of rules (algo- 
rithms) to compute probability distributions [5,8,9]. Therefore, we may wonder what kind 
of algorithm(s) we need to perform an event-based simulation of the experiments mentioned 
above without using wave functions. Evidently, this formulation rules out any method based 
on the solution of the (time-dependent) Schrodinger equation. We have to step outside the 
framework that quantum theory provides. 

We formulate the simulation of quantum processes in terms of events, messages, and units 
that process these events and messages. In terms of the experiments of Grangier et al. [1], the 
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photon carries the message (a phase), an event is the arrival of a photon at one of the input 
ports of a beam spHtter, and the beam sphtters are the processing units. The essential feature 
of a processing unit is its ability to learn from the events it processes. We use a standard 
linear adaptive filter [10] for this purpose. The processing unit sends a message (carried by 
a photon) through an output port that is selected randomly according to a distribution that 
is determined by the current state of the processing unit. A processing unit that operates 
according to this principle will be referred to as a stochastic learning machine (SLM). Here 
the term stochastic does not refer to the learning process but to the method that is used to 
select the output channel that will carry the outgoing message. The learning process itself is 
deterministic. Elsewhere we describe processing units that are fully deterministic and are not 
based on linear adaptive filters, but can be used for exactly the same purpose and also for 
simulating many-body quantum phenomena [11-14]. 

By connecting an output channel to the input channel of another SLM we can build 
networks of SLMs. As the input of a network receives an event, the corresponding message 
is routed through the network while it is being processed. At any given time during the 
processing, there is only one input-output connection in the network that is actually carrying 
a message. The SLMs process the messages in a sequential manner and communicate with 
each other by message passing. There is no other form of communication between different 
SLMs. Although networks of SLMs can be viewed as networks that are capable of unsupervised 
learning, they have little in common with neural networks [15]. 

Beam splitter. ~ According to quantum theory [9], the amplitudes (6o, bi) of the photons 
in the output modes and 1 of a beam splitter (see Fig. ^ are given by 

(bo\^J_fao + iai\^J_/'li\/'ao\ 

where the presence of photons in the input modes or 1 is represented by the amplitudes 
(ao,ai) [1,16-18]. Here we construct a SLM that acts as a beam splitter, not by calculating 
the amplitudes according to Eq. but by processing individual events (photons). First we 
describe the various components of the machine. Then we demonstrate that it acts as a beam 
splitter. 

A schematic diagram of the SLM is shown in Fig.^ We label events by a subscript n > 0. 
At the {n + l)-th event, the SLM receives a message on either input channel or 1, never 
on both channels simultaneously. If the event occurs on channel (1) we call it an event 
of type (1). Every message is a two-dimensional unit vector y„+i — (yo,n-i-i, ?/i,n+i)- This 
message represents the time-of- flight, that is the length of the optical path of the photon when 
it travels from source to beam splitter, from beam splitter to beam splitter and the like. In 
quantum theory, this information is encoded in the phases of the complex numbers ap and oi. 

The first stage (see Fig. of the SLM (called front-end) stores the message yn-i-i in its 
internal register = (yo^k,Yi.k) where k — (1) if the event occurred on channel (1). 
The front-end also has an internal two-dimensional vector x — {xo,xi) with the additional 
constraints that a;^ > for i = 0, 1 and that xq+xi = 1. Here and in Fig.Q] we have omitted 
the event-label n from the internal register and internal vector x to simplify the notation. 

After receiving the (n + l)-th event on input channel fc = 0, 1 the internal vector is updated 
according to the rule 

Xi,n+i = aa;i,ri + (1 - a)(5i,fc, (2) 

where < a < 1 is a parameter that controls the learning process and is discussed later. By 
construction Xi^n+i > for i = 0,1 and a;o,„+i + Xi^n+i = 1- Hence, the update rule Eq. ^ 
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Fig. 1 - Diagram of a SLM that performs an event-by-event simulation of a single-photon beam 
splitter (BS). The solid lines represent the input and output channels of the BS. Dashed lines indicate 
the flow of data within the BS. 



preserves the constraints on the internal vector. These constraints are necessary if we want to 
interpret the Xk,n as (an estimate of) the probability for the occurrence of an event of type k. 
The solution of Eq. (0) reads 

n-l 

x„ - a"xo + (1 - a) ^ V,+i , (3) 

i=0 

where x„ = (so.n, and Xq denotes the initial value of the internal vector. The input 

events are represented by the vectors v„+i = (1,0)-'" or v„+i (0, 1)"^ if the n + 1-th event 
occurred on channel or 1, respectively. Let pq {pi = 1 — po) be the probability for generating 
statistically independent events of type (1) and let (x)„ — (xi + . . .+x„)/n denote the Cesaro 
mean of the sequence {xi, . . . , x„}. Then, if < a < 1, we have lim„^oo(x)n = {po, ^— Po)^ ■ 
Therefore, the front-end "learns" the probabilities for events and 1 by processing these 
events in a sequential manner. 

The second stage of the SLM in Fig. Stakes as input the values stored in the registers Yq, 
Yi, X and transforms this data according to the rule 
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Fig. 2 - Simulation results for the beam splitter shown in Fig. Q Input channel receives 
(cos i/)o, sin ^o) with probability po- Input channel 1 receives (cos sin-!/)!) with probability pi = 
1 — po- Each data point represents 10000 events (A'o + A''i = 10000). After each set of 10000 events, 
two uniform random numbers in the range [0, 360] are used to choose the angles tpo and ^pi. IVIarkers 
give the simulation results for the normalized intensity No/{No + Ni) in output channel as a function 
of = t/'o — tpi. Open circles: po — 1; Bullets; po ~ 0.5; Open squares: po = 0.25. Lines represent 
the results of quantum theory (see Eq. 0)- Left: a = 0.98; Right: a = 0.25. 



where we have omitted the event label (n + 1) to simplify the notation. Using two complex 
numbers instead of four real numbers Eq. 10} can also be written as 



(5) 



Identifying uq with lo,o\/^ + iYi.o^/xo and ai with Yq^i^^ + iY\^\yfx\ it is clear that the 
transformation Eq. (|SJ) plays the role of the matrix- vector multiplication in Eq. iQJ . 

The third stage (see Fig. ^ of the SLM (called back-end) responds to the input event 
by sending a message w„-|_i = (lo,o\/a^ ^ Y\^\^fx{^ Yi^^\^fx\ -I- Yx^^/xq)! \f2 through output 
channel if + w\ < r where < r < 1 is a uniform random number. Otherwise, 

since w^^^ +z,^j_,_^ = 1, the back-end sends the message z„+i — iyQ,\^/x\ — Y\^Q^fxQ, Yq^q^/x^^ 
Y\^\sjx\) I \f2 through output channel 1. 

In Fig. 13 we present results of event-based simulations using the SLM depicted in Fig.^ 
Before the simulation starts we set xq = (xq.o, 2;i,o) = 1 — 'r) where r is a uniform random 
number in the interval [0,1]. We also use uniform random numbers to generate two two- 
dimensional unit vectors to initialize the registers Yq and Yi . The simulation results of the left 
and right panel were obtained for a — 0.98 and a — 0.25, respectively. Input channel receives 
(cos "00 1 sin ^/jq) with probability po- Input channel 1 receives (cos ^/^i , sin ) with probability 
Pi = 1 — Po- III the quantum-theory setting, this corresponds to the input amplitudes ag = 
^J^e^'^^ and = — poe^^^ . The data in Fig. |21 show the normalized intensity Na/{NQ + 
Ni), where Nq (Ni) denotes the number of events of type (1), as a function oi 4> = ip^ — ipi. 
Each data point in Fig. [21 represents a simulation of 10000 events (A'o -I- A^i = 10000). For 
each set of 10000 events, two uniform random numbers in the range [0, 360] determine the two 
angles ■00 and 0i • From Fig. 12 it is clear that the SLM-based beam splitter reproduces the 
probability distributions 

i, i2 1 + 2v/po(l - Po) sin(0o - -01) |, |2 1 - 2 Vpo(1 - Po) sin(V'o - ipi) 

M = ^ , \bi\ = , (6) 
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Fig. 3 - Diagram of a SLM network that simulates a single-photon Mach-Zehnder interferometer 
on an event-by-event basis [14]. The SLM network consists of two BS devices (see Fig. ^ and 
two passive devices R{(j)o) and R((j)i) that perform plane rotations by cj>o and <j)\, respectively. The 
number of events Ni in channel i = 0, . . . , 3 corresponds to the probability for finding a photon on 
the corresponding arm of the interferometer. 



as obtained from Eq.QJ. 

Role of the control parameter a. - From Eq. (O, we see that a controls the speed of 
learning. Furthermore, it is evident that the difference between a constant input to a SLM 
and the learned value of its internal variable cannot be smaller than 1 — a. In other words, a 
also limits the precision with which the internal variable can represent a sequence of constant 
input values. On the other hand, the number of events N has to balance the rate at which the 
SLM can forget a learned input value. The smaller 1 — a is, the larger the number of events 
has to be for the SLM to adapt to changes in the input data. Simulations (results not shown) 
show that the error decreases as l/y/N, as expected from elementary statistics. 

In general, for a fixed number of events, decreasing a leads to an increase of systematic 
errors. For instance, if a = 1/4 the maximum normalized intensity (at </> = 90) in output 
channel is about 0.8, as shown in Fig. El (right). This is easy to understand: In our event- 
by-event approach, interference is the result of learning by the SLMs. If a is close to one, 
the SLM learns slowly but accurately. If a decreases, the SLM can adapt faster to changes 
in the input data but it also forgets faster. For very small a, x„ is always close to (1,0)'^ 
or (0, 1)"^ and it becomes impossible to mimic quantum interference. Thus, a determines the 
visibility [1] of the interference effects. In our simulation approach, interference is the result 
of learning (which requires some form of memory) . 

Mach-Zehnder interferometer. - In quantum physics [9] , single-photon experiments with 
one beam splitter provide direct evidence for the particle-like behavior of photons [1,5]. The 
wave mechanical character appears when one performs single-particle interference experi- 
ments. We now describe a SLM network that displays the same interference patterns as those 
observed in single-photon Mach-Zehnder interferometer experiments [1]. This also proves that 
the messages generated by the SLMs preserve the phase information that is essential if the 
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Fig. 4 - Simulation results for the SLM-network shown in Fig. |2| Input channel receives 
(cos V'o, sin ^o) with probability one. A uniform random number in the range [0,360] is used to 
choose the angle i/jq. Input channel 1 receives no events. Each data point represents 10000 events 
(A^o + Ni = N2 + N3 ^ 10000). Initially the rotation angle (po = and after each set of 10000 events, 
(po is increased by 10° . Markers give the simulation results for the normalized intensities as a function 
of 4> = 4>o ^ 4>i ■ Open squares: No/{No + Ni); Solid squares: ^2/(^2 + N3) for = 0; Open circles: 
A2/(iV2 + N3) for (^1 = 30°; Bullets: ^"2/(^2 + A3) for </)i = 240°; Asterisks: A3/(A2 + A3) for 
01 = 0; Solid triangles: N-j,/{N2 + A3) for 01 = 300°. Lines represent the results of quantum theory 
(see Eq. ©). Left: a = 0.98; Right: a — 0.25, see the section on the role of a. 



system is to exhibit quantum mechanical behavior. 

The schematic layout of the SLM network is shown in Fig. El Not surprisingly, it is exactly 
the same as that of a real Mach-Zehnder interferometer. The SLAI network described before is 
used for the two beam splitters in Fig. |3| The phase shift is taken care of by a passive device 
that performs a plane rotation. Clearly there is a one-to-one mapping from each physical 
component in the interferometer to a processing unit in the SLM network. The processing 
units in the SLM network only communicate with each other through the messages (photons) 
that propagate through the network. 

According to quantum theory [9], the amplitudes (69,^1) of the photons in the output 
modes {N2) and 1 {N^) of the Mach-Zehnder interferometer are given by [1,16-18] 




Note that in a quantum mechanical setting it is impossible to simultaneously measure {Nq / {No+ 
Ni), Ni/{No + Ni)) and {N2/iNo + Ni), Ns/iNo + Ni)): Photon detectors operate by ab- 
sorbing photons. In our event-based simulation there is no such problem. 

In Fig. 0] we present a few typical simulation results for the Mach-Zehnder interferometer 
built from SLMs. We assume that input channel receives (cos ipo, sin ipo) with probability one 
and that input channel 1 receives no events. This corresponds to (ag : ) = (cos ^/jq+i sin V'o , 0) . 
We use a uniform random number to determine ipQ. The data points are the simulation results 
for the normalized intensity Ni/{NQ + Ni) for i=0,2,3 as a function of = 00 — 0i. Lines 
represent the corresponding results of quantum theory [9]. From Fig. ^ (left) it is clear that 
the event-based processing by the SLM network reproduces the probability distribution 

|6o| =sm( ), [foil =cos( ), (8) 
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as obtained from Eq.([7I). The results of Fig. 0] (right) for a = 1/4 corroborate the analysis of 
the role of a, presented earlier. 

Discussion. - We have proposed a new procedure to construct algorithms that can be 
used to simulate quantum processes without solving the Schrodinger equation. We have shown 
that single-particle quantum interference can be simulated on an event-by-event basis using 
local and causal processes, without using concepts such as wave fields or particle- wave duality. 
Our results suggest that we may have discovered a procedure to simulate quantum phenomena 
using causal, local, and event-based processes. There is a one-to-one correspondence between 
the parts of the processing units and network and the physical parts of the experimental 
setup. Only simple geometry is used to construct the simulation algorithm. In this sense, the 
simulation approach we propose satisfies Einstein's criteria of realism and causality [5] . Our 
approach is not an extension of quantum theory in any sense and is not a proposal for another 
interpretation of quantum mechanics. The probability distributions of quantum theory are 
generated by a particle- like, causal learning process, and not vice versa [7]. 
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